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Abstract. We present results from our numerical simulations of collapsing 
massive molecular cloud cores. These numerical calculations show that massive 
stars assemble quickly with mass accretion rates exceeding 10 -3 MQyr -1 and 
confirm that the mass accretion during the collapsing phase is much more effi- 
i <~| cient than predicted by selfsimilar collapse solutions, M ~ c 3 /G. We find that 

Oh' during protostellar assembly out of a non-turbulent core, the mass accretion 

reaches 20 — 100c 3 /G. Furthermore, we explore the self-consistent structure of 
bipolar outflows that are produced in our three dimensional magnetized collapse 
^ ■ simulations. These outflows produce cavities out of which radiation pressure can 

| be released, thereby reducing the limitations on the final mass of massive stars 

formed by gravitational collapse. 

Additional enhancement of the mass accretion rate comes from accretion 
along filaments that are built up by supersonic turbulent motions. Our numerical 
calculations of collapsing turbulent cores result in mass accretion rates as high 
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Our understanding of how massive stars form, while still far from complete 
has recently made significant strid es (see also Harold Yorke 's contribution in 
this proceedings and the review by IZinnecker Yorkd 120071 b Low mass stars 



accrete the bulk of the ir mass through their circumstellar disks before nuclear 



burning turns on (e.g. JShu et"aDll987l b Massive stars on the other hand, have 



Kelvin-Helmholtz time scales that are smaller than the dynamical time so that 
young massive stars s tart to burn their nuclear fuels while still accreting the 
surrounding gas (e.g., Yorkd 2002. 2004, and references herein). Infall, or flow 



of gas through a surrounding disk, therefore faces a major obstacle in the form 
of the radiative pressure that such massive stars will produce as they are still 
forming. 

Early spherical accretion models suggest that the resulting radiation pres- 
sure could limit the final mass of the star to ~ 40 if the accretion ra te is not 
high enough ( > 10~ 3 M Q yr" 1 ) (lKahnl ll974l: IWolfire & Cassinelli| Il987lh More 



recent two dimensional simulations bv lYorke k, Sonnhalterl ([2002) showed that 
the limitations on the final mass of the massive star can be relaxed if accretion 
through the protostellar disk is included in models of massive star formation 



Still, even this particular result yields mass limits of ~ 43 Mq. iKrumholz et al 
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Figure 1. Shows the time evolution of the column density and mass ac- 
cretion of a collapsing magnetized core. We assume the beginning of the 
collapsing phase at t = to when the column density reaches Score ~ lgcm~ 2 
and the pressure is P/ks ~ 10 s K cm -3 . At t = to the mass density in 
the core is pcorc — 2.44 x 10 ^gem 3 corresponding to a free fall time of 
iff = 1.35 x 10 4 yrs. At t = to + 1.1% the core density reached already 
10 4 g cm~ 2 . The initial profile o f our simulation is marked by t = 0. [Adapted 
from lBaneriee fc Pudritd (|2007j )] 



( 20051 ) argued that the effects of radiation pressure are limited by the escape of 
radiation through an outflow cavity. These results were based on simulations, 
using Monte-Carlo-diffusion radiative transfer models where the outflow cavities 
are parameterized by varying opening angles. 

Alternative ly, massive stars co uld form through coalescence of intermedi- 



ate mass stars ( Bonnell et al. 19981 ). This formation process would be starkly 



different from the formation of low mass stars which assemble quickly through 
accretion of the molecular gas. So far, observations of the intermediate state of 
massive star formation are rare and difficult to obtain. Nevertheless, a few mas- 
sive objects which show evidence for an ongoing accretion process are known b y 
now (IChini et alJl2004 : iPatel et al.ll2005l : IChini et al.ll2006l : iBeltran et alJl200d ). 

Here we summarize our recent effort studying the collapse of massive cloud 
cores which final state will be one or more massive stars. These cores are either 
modeled as supercritical hydrostatic spheres or taken from simulations of super- 
sonic turbulence in which the unstable cores are formed in shock compressions. 
The numerical setups and the technical details which also inclu des our descrip- 
tion of var ious cooling and heating p rocesses can be found in IBaneriee et al 
(j2006T ) and IBaneriee fc Pudritz! <j2007T ) . 



We find that Bonnor-Ebert models for collapsing magnetized cores result 
in enormous accretion rates and outflows that are driven by toroidal magnetic 
pressure that is built up in the massive disk. When filamentary structure is 
added to this picture, infall becomes even stronger as gas is funneled from large 
scales into the filament, and then down onto the forming disk. 



2. High accretion rates during the collapse phase 

The collapse of spherical cloud cores h as been studi e d by many au t hors , both 



analytically and numerically (see, e.g., Larson 196S ; Penstonl Il969j; IShu 19771: 



Hunterlll977l : IWhitworth fc Summersll985l : lFoster fc Chevalierlll993l : lHennebelle etal 
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120031 : iBaneriee et alJl2003 ). In the case of an initial singular isothermal sphere 
(SIS) there exist an elegant self-similar solution to this problem where an ex- 
pansion wave travels from the (singula r) center outward with the speed of 
sound initiating an ins ide-out collapse ( Shul Il977l ). But, as pointed out by 
Whitworth et all ( 19961 ) a SIS configuration is unnatural be cause - among other 
difficulties - its collapse can not produce binaries (see also |PringleHl989r i, The 



collapse of non-singular cores proceeds differently than singular spheres and has 
distinguishab l e implications : First - as demonstrated by many authors (e.g. , 
Larson 1969; iPenstonl Il969l : iFoster fc Chevalier! Il993l : iHennebelle et al.1 [2003; 
Baneriee et alJl2004T ) - the collapse proceeds from outside-in rather than from 
inside-out and the density maintains a flat profile at the core center, where 
the core size is of the order of the local Jeans length at every epoch (see also 
Whitworth &: Summer i l!985l . for a summary of analytic solutions). Second, the 



radial distribution of the infall velocity peaks at the edge of the flat density cor e 
and falls off quickly towards the center (e.g., see Fig. 3 in IBaneriee etaD 120041 ). 
The velocity becomes also supersonic as the core density increases and the size of 
the flat region shrinks. Recent observations of a pre-Class object show that the 
collapse proceeds supersonically, strongly supporting a Larson- Periston- t ype col- 
lapse rather than an expansion wave- type collapse ( Furuva et al.l 120061 ) . Third 
- and of great importance - the mass accretion in the non-selfsimilar case is 
much higher than predicted from the selfsimilar collapse of a singular isother- 
m al sphere. We find that th e mass accretion in the early phase of the collapse 
is (|Baneriee Pudritd [20071 ) 



M w 20 - 100 — 



(1) 



(c is the isothermal sound speed and G is Newton's gravitational constant). Note 
the selfsimilar SIS collapse gives a mass accretion of only 0.96 c 3 /G. Typical 
values of the sound speed in cold cloud cores (T ~ 20 K) are of the order of 
a few 10 m sec - which gives c 3 /G ~ 10~ 6 - lO~ 5 M0yr -1 . Our resu lt from 
numer ical simulations are in agree ment with the early analytic results of iLarsonl 
(119691 s ) and IPenstonl ([19691 ) (see also lHuntedll977l : IWhitworth fc Summerslll985l ). 

The remarkable point is that the high accretion rates are achieved even 
without initial turbulence and during the isothermal phase of the collapse. The 
main reason for the high accretion rate in this idealized case is the supersonic 
infall velocity, v r , close to the peak density. Even moderate Mach numbers, A4, 
of 2 — 3 (the Larson asymptotic Mach number is A4 ~ 3) enhances the mass 
accretion rates relative to c 3 /G because 



VG -> vi/G = M 6 c 6 /G 



(2) 



in the supersonic limit. Additionally, the core is continually embedded in a high 
pressure environment. 

Studies of the long t erm evolution (i.e., beyon d one dynamical time) of a 
collapsing BE-Sphere by IFoster Sz Chevalier] ( 19931 ) showed that the mass ac- 
cretion is not constant, but decreasing with time after it reached a peak value. 
Beyond this turnover point the envelope is drained of gas resulting in a decreas- 
ing accretion rate. This situation might be different for the collapse of more 
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Figure 2. Comparison of the radial infall velocity for the three differ- 
ent cases, hydromagnetic (mag), pure hydro (hydro), and pure isother- 
mal (iso) simulation at the time they reached the same core density £ = 
5 x 10 3 gcm -2 . The left panel shows the infall velocity in km/sec and com- 
pared to the initial isothermal sound speed, q so , and the right panel shows 
the s ame quantity measured w ith the local sound speed, q 0C al- [Adapted 
from lBanerjee & Pudritel (|2007l )] 



realistic cloud cores which are not isolated objects but are surrounded by a 
clumpy medium whose accretion might sustain high accretion rates for a longer 
time. 

At the time when the column density reaches £ ~ 1 g cm -2 the mass accre- 
tion rate becomes ~ 10 -4 Mq yr _1 and the surrounding pressure is 10 s Kcm~ 3 fcR. 
These result s are in agreement with the turbulent core collapse model of lMcKee Sz Tan 
( 20o3. 120031 ) who showed that this high pressurized, compact cloud cores accrete 
gas with accretion rates as high as 10 -3 Mq yr _1 . We find that accretion rates of 
this order (~ 10~ 3 Mq yr _1 ) are reached within only ~ 14.7 x 10 3 yrs which cor- 
responds to 1.1 dynamical times where core column density reaches ~ 10 4 gcm~ 2 
(see Fig. [I]). As long as the core stays isothermal (the efficient cooling regime) 
the pressure scales with the column density as P oc £ 2 (i.e., the pressure profile 
is close to R~ 2A ) and slightly steeper in the inefficient cooling regime where the 
temperature rises during the collapse. The continuous high external pressure 
and the supersonic infall velocity maintains the high mass accretion. 

In Fig. [T] we also show the mass accretion scaled to the quantity cf so /G (ci so 
is the initial isothermal sound speed) from which one can see that accretion is 
much more efficient (by a factor of ~ 100) than expected from the collapse of a 
SIS. Again, the reason for the high accretion rates is the supersonic infall of gas. 
Fig. [2] shows the infall velocities for a magnetized (mag), pure hydro dynamical 
case with cooling (hydro), and an isothermal collapse. Here, the infall velocities 
reach up to 1.5 km see -1 which corresponds to four times the initial sound speed. 
In all cases the infall velocities also exceed the local sound speed where the local 
Mach numbers vary between 2 and 2.5. 

These results show that high mass accretion rates are a natural result of 
the early collapse phase of (non-singular) collapsing cloud cores with flat-topped 
density profiles. Turbulent driving for the rapid assembly of massive stars is not 
a necessary ingredi ent, but it does further enh ance the accretion rate as shown 
in analytic models ( McKee Tan 20021. 2003) and ev en more strongly in fully 



realized 3D turbulence (see, e.g., Banerjee et al.|[2006l ). 
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Figure 3. Close-up 2D snapshots of the disk region (perpendicular to the 
disk in the yz-plane) in the magnetized case. The left panel shows the sit- 
uation at t = 1.45 x 10 4 yrs (1.08%) into the collapse and before the flow 
reversal and right panel shows the c onfiguration 188 yrs later w hen the out- 
flow is clearly visible. [Adapted from iBaneriee fc Pudritd (|2007l ).] 



3. Early outflows during massive star formation 

Very little is known observationally about the influence of outflows on the early 
assembly of massive stars. On the one hand hand, they could reduce the mass 
accretion onto the massive (proto)star if the outflow carries a substantial mass. 
On the other hand, early outflows provide a natural anisotropy of the accret- 
ing gas which results in low density cavities. Such cavities are like funnels 
out of which the radiation from the already active star can escape. With- 
out such a pressure-re lease valve, trapped radiative flux would halt the infall. 
Krumholz et al. ( 20051 ) studied the effect of outflow- funnels using a Monte Carlo 



radiative-transfer method which shows that radiation pressure is greatly reduced 
by radiation escaping the outflow cavities. 

Magnetic fields coupled to the protostellar disk can be the driving power for 
such outflows. A variety of self-consistent simulations of collapsing magnetized 
cloud cores show that outflows are la unched if the toroidal magnetic field pres- 
sure o vercomes the gravita t ional force ( Tomisakall 19981 . 120021 ; iMatsumoto Tomisaka 



20041 : iMachida et alJl2003 : IBaneriee k Pudritzll2006l ) . Such early-type outflow s 



can be understood in terms of a growing magnetic tower ( Lvnden-Belll 120031 ) : 
The rotating (proto)disk generates a strong toroidal field component by winding 
up the threading field lines. The resulting magnetic pressure is in local equilib- 
rium with the gravitational force and the ram pressure of the infalling material. 
But every new rotation increases the toroidal field component thereby shifting 
the equilibrium location (which is characterized by a shock) towards higher lat- 
itudes. The result is an inflating magnetic bubble in which material is lifted off 
the disk. 

We explicitly demonstrate the launch of the outflow in the 2D cuts through 
our 3D data, shown as the images in Fig. [3] which show the collapse state at 
different times. These close-ups (~ 20 AU) of the early stage of the outflow 
launching show the collapsing stage shortly before the flow reversal (left panel) 
and the onset of the outflow (less than 200 yrs later, right panel). 
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Figure 4. Collapse of a elongated tur ubulent core which is b uilt up in a 
supersonic turbulent molecular cloud (see lTillev fc Pu dritz 2005). Shown are 
the density isosurfaces at different scales: the filamentary structure (2400 AU), 
the protostellar disk within t he filament (150 AU) , a nd the inner region of the 
disk (37 AU). [Adapted from lBaneriee et all (|2006t ).] 



4. Massive star assembly in supersonic turbulence 



Our current understanding of star formation is that stars are a natural conse- 
quence of supe rsonic turbulence within s elf-gravitating, molecular clouds (see re- 

cent r eviews by McKee k, Ostrikerll2007l ; iBallesteros-Paredes et al.ll2007l ; lElmegreen k, Scaio 
20041 : iMac Low fc Klessenll2004l ). 

Supersonic turbulence is observed in most if not all giant molecular clouds 
(GMCs) and is important because it rapidly sweeps up large volumes of gas and 
compresses it into systems of dense filaments. Observational surveys confirm 
that filamentary substructure characterizes the internal organization of molec- 
ular clouds on many scales. It is c learly evident in the Orion A molecular 
cloud (e.g., iJohnstone Bally 1999) where, in addition to showing the obvi- 
ous integral-filament shaped structure, one also sees smaller structures o f 1.3 pc 
in scale. JCMT submillimeter studies of Orion B ( Mitchell et al. 200 ll ) reveal 
a plethora of filamentary structures an d their embedded cores. This pattern is 
also seen in the observational studies bv lFiege et al.l (2004), wherein a few bright 
cores are seen to be embedded in a larger, isolated filamentary structure. Sim- 
ilar results are also seen in more embedded regions such as the Lupus 3 cloud 
where strong links b etween filaments and emerging star clusters are observed 
(|Teixeira etal.1 120051 ) 

Following the collapse of a massive elongated core which is for med within a 
super sonic environment (the initial data where taken from one of the lTillev &i Pudritz 
(2005) simulations), we explore the multi-scale nature of the formation of fila- 
mentary large scale structure and how this controls the formation of protostellar 
disks and stars in turbulent cluster-forming clumps. We follow the filament that 
forms the first massive star in our simulation, and find that it has high density 
and dynamic pressure that is maintained by continued inflow and shocking into 
the filaments. The collapse of this dense material along the filament and into 
the disk and star is the key to understanding this problem. Lower mass stars 
in this picture also form in filamented structure, but which are less compressed 
and on smaller scales. Our simulations trace how the most massive star-forming 



Filaments, Collapse <fe Outflows 7 

region is assembled and how collapse and the formation of protostellar disks can 
occur, by resolving the local Jeans length down to scales approaching that of 
the protostar itself. 

We find that filaments play a dominant role in controlling the physics, accre- 
tion rate, and angular momentum of the much smaller-scale accret ion disk that 
forms within such collapsing structures (see also lBalsara et al.ll200ll ). Large-scale 
filamentary flows sustain accretion rates that are orders of magnitude greater 
than both (i) the naive scalings that are derived from the virial theorem ap- 
plied to uniform, 3D media, or even (ii) the collapse of isolated Bonner-Ebert 
spheres. The large-scale structuring of molecular clouds into filaments therefore 
has profound effects on the rate of formation of disks and stars. 

In Figure U] we show the 3D structure of the collapsing filament at different 
scales which summarizes our results. The large-scale filament is about 2600 AU 
in length and one sees the converging accretion flow on the 150 AU scale view 
presented in the middle panel. On this intermediate, 150 AU scale, the disk 
grows within a filament that is developing out of a sheet-like structure. The 
transfer of material from the sheet into the filament is associated with net an- 
gular momentum - and tends to resemble the growth of a large-scale, spinning 
vortex. The disk, forming on yet smaller scales, acquires this angular momen- 
tum from this larger scale process. We see the actual disk at this time in the 
right panel, on a scale of only 37 AU in diameter. 

One of the major results of this the filamentary collapse is the very high 
accretion rates (M ~ 10 _2 MQyr _1 ) due to the supersonic gas flow onto the 
protostellar disks. These rates are 10 3 times larger than predicted by the collapse 
of singular isothermal spheres and exceed the accreti on rates necessary to squeez e 
the radiation field of the newly born massive star (jWolfire Cassinellil \l98l\ ) . 



We find that a reasonable scaling for filamentary accretion is M ~ / c 3 /G where 
the prefactor should scale as / ~ Ai 3 , i.e., that one replaces the sound speed 
in the virial formula with the free-fall speed. Even this does not quite account 
for the full effect, and there is an additional factor of order a few that likely 
is geometric in nature and accounts for filamentary, rather than spherical infall 
geometry. Filamentary structure, therefore, plays a very significant role in star 
formation by funneling flows, gathered from large sales, along the filaments and 
into the disks. 

Our future work will pursue these insights by including the diffuse radiation 
from the accretion process and ionization feedback from the massive star. This 
will show whether our results are relevant for more realistic environments. 
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